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Abstract. The Time-Dependent Orienteering Problem (TDOP) is a 
generalization of the Orienteering Problem where graph weights vary 
with time. It has many real life applications particularly associated with 
transport networks, in which travel time between two points depends on 
the moment of start. The paper presents an evolutionary algorithm with 
embedded local search operators and heuristic crossover, which solves 
TDOP. The algorithm was tested on TDOP benchmark instances and in 
most cases achieved optimal or near optimal results clearly outperform¬ 
ing other published methods. 

Keywords: Time-Dependent Orienteering Problem, evolutionary algo¬ 
rithm, metaheuristic, local search operators 


1 Introduction 

The Orienteering Problem (OP) can be formulated on a graph in which each 
vertex has assigned a profit and each edge has a cost. The goal of the OP is 
to find a path between given points limited by total cost which maximizes to¬ 
tal proht of visited vertices. The OP has many practical applications including 
tourist trip planning [17], transport logistics and others. For example in tourist 
trip planning each attraction (point of interest - POI) has some profit assigned. 
In this case the purpose is to find as attractive trip as possible with limited 
duration. It can by solved by OP algorithms. 

In this paper the Time-Dependent Orienteering Problem (TDOP) [15] is tack¬ 
led. TDOP is a generalization of OP in which graph weights (indentified with 
travel times) vary with time. In practical situations, usually travel time between 
two places depends on the moment of travel start. Examples of time-dependent 
networks are public transport networks (travel times determined by timetables) 
as well as other transport networks (travel times determined by time of a day 
and traffic). Determining most profitable paths in such networks (i.e. most at¬ 
tractive trips for tourists) comes down to solving TDOP. 

To solve the TDOP author proposed an evolutionary algorithm with local 
search methods embedded. It was based on the algorithm which author devel¬ 
oped for classic OP obtaining solutions close to optimal for benchmarks up to 
a few hundred vertices [27]. Some modihcations were applied to the algorithm: 
allowance of infeasible solutions (too long paths) in a population and adaptation 
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of operators to time dependency. 

The article is organized as follows. Literature review is presented in section 2. 
Detailed mathematical problem definition is given in section 3. In section 4 there 
is a description of proposed evolutionary algorithm. Section 5 contains compu¬ 
tational results of the algorithm and comparisons with other metaheuristics. 
Finally, in section 6 conclusion and further research are included. 


2 Literature review 

The Orienteering Problem was introduced by Tsiligirides [1]. The OP (and its 
generalizations) is an NP-hard optimization problem [2] and exact algorithms for 
OP can be very time-consuming for larger graphs. The exact algorithms usually 
applied for the OP are the branch-and-cut methods [7] [5]. However, because 
of its computational difficulty, most OP papers are devoted to heuristics and 
metaheuristics. Various approaches were based i.a. on greedy and randomized 
construction of solutions [1] [2] [20], local search methods [4] [12], tabu search [6], 
ant-colony optimization [13] and genetic algorithms [8]. Author’s previous OP re¬ 
search concentrated on development of evolutionary algorithms with local search 
operators, which proved successful on benchmark instances [28] [27] [25] [21]. 

The Time-Dependent Orienteering Problem (TDOP) is relatively new and 
most of papers regarding TDOP were published in current decade. Most of 
them emphasize practical application of TDOP especially in transport networks. 
Fomin et al. [9] presented first approach to the problem - they proposed (2 -|- e)- 
approximation algorithm for TDOP. The problem version tackled in the article 
was simplified because each node had the same profit (equal to 1). LI [15] de¬ 
veloped an exact method basing on dynamic programming and solved small 
instances of TDOP. He used discretized time and there were 1440 time states 
(number of minutes in a day). It was the first article about classic TDOP (nodes 
profits can be arbitrary). 

The TDOP solutions can be practically used in networks of tourist attrac¬ 
tions connected by means of public transport. A few papers concentrated on 
TDOP application in such networks. In paper [18] first such an application was 
presented - algorithm was tested on POI and public transport network of Tehran. 
The authors used a different name to the problem and routes utility was com¬ 
puted in a different way to classic OP but the problem was closely related to 
TDOP. 

Garcia et al [19] for the first time applied classic TDOP version in a public 
transport network - experiments were conducted on POI and public transport 
network of San Sebastian. Their solution also took into account time-windows 
and multiple paths (Time-Dependent Team Orienteering Problem with Time 
Windows). In order to solve the problem they used Iterated Local Search method 
(ILS). However, they used some simplifications: computations performed on av¬ 
erage daily travel time (first approach) and assumption of periodic public trans¬ 
port timetables (second approach). 

Recently one of the first approaches which uses real time-dependent weights 
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in public transport network was presented [26]. Experiments were conducted on 
POI and public transport network of Athens. Authors proposed 20 topologies 
and 100 tourist preferences combining into 2000 different test cases. Authors 
introduced two fast heuristics (TDCSCR and TDSICSCR), which based on it¬ 
erated local search and vertex clustering, and made comparisons with ILS. 

Another solution of TDOP tested on real travel times of POI and public 
transport network (city of Bialystok) was proposed by the author of this paper 
[24]. The algorithm based on local search procedures and was tested in three vari¬ 
ants which processed network weights in different ways. Time-dependent variant 
performed computation on real travel times, mean variant operated on static 
network of daily averaged travel times (like in approaches described above) and 
hybrid variant used information from both approaches. 

Verbeeck at al [22] developed new benchmark instances for TDOP. These 
networks model street traffic. The authors proposed ant-colony approach, tested 
it on mentioned benchmarks and compared with optimal solutions (obtained by 
them). The tested metaheuristic achieved high quality results in short execu¬ 
tion time. Gunawan et al [23] modified Verbeeck’s benchmarks (discretization 
of time) and tested a few approaches: greedy construction, iterated local search 
(ILS), adaptive ILS and variable neighborhood descent. The best solutions were 
generated by adaptive ILS metaheuristic. 


3 Definition of the Time-Dependent Orienteering 
Problem 


Let G = (R, E) be a directed, weighted graph. In this paper each weight between 
vertices i and j (i,j G V) is identified with travel time between these vertices 
and is a function Wij{t) dependent on the moment of travel start t. Each vertex 
i has a nonnegative profit pi and a nonnegative visit time r^. Given the time- 
dependent graph, moment of start to, time limit T^ax, start and end vertices (s 
and e) the purpose of TDOP is to find a path from s to e starting at time to 
which maximizes total profit of visited vertices and its total cost (travel time) 
does not exceed T^ax- TDOP can be formulated as Mixed Integer Programming 
problem. Let Xij is 1 iff a path contains direct travel from i to j and 0 otherwise. 
Let Pi is 1 iff a path contains vertex i and 0 otherwise. Let ti be a time of arrival 
at vertex i - this function is defined only for vertices contained in a path. The 
purpose of TDOP is to maximize formula 1 without violating constraints given 
by formulas 2-8. Eormula 2 limits total time budget. Eormula 3 indicates that the 
path starts in s and ends in e. Constraint 4 means that each vertex can be visited 
only once while formula 5 is a relation between variables x and y. Eormulas 6 and 
7 are used to compute arrival times at visited vertices. Constraint 8 guarantees 
that there are no sub cycles. 
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4 Proposed algorithm description 

To solve the TDOP author proposed an evolutionary algorithm with local search 
methods embedded (during mutation), disturb operator, 2-point crossover and 
deterministic crowding as selection mechanism. Operators have two forms (ran¬ 
dom and heuristic) - their frequency and specificity are controlled by algorithm 
parameters. Infeasible solutions (too long paths) can be included in a population 
as well and they are penalized by fitness function. A path representation is used 
in the proposed algorithm - each gene in a chromosome is associated with visited 
vertex in a corresponding path. After random initialization of paths, procedures 
of crossover, selection, mutation and disturbance and applied repeatedly. Evo¬ 
lution terminates after Ng generations or earlier if there is no improvement in 
the last Cg generations. The result is the best feasible path obtained during 
algorithm run. 

4.1 Evalutation 

Let S' be a path with total profit ps and total cost (travel time) eg. Fitness 
function of a feasible path (not exceeding Tmax) is equal to its profit. If path S 
is too long its fitness is: fitness{S) = ps ■ ( ^^^^ )^. Fitness function decreases 
as solution is farther from feasibility boarder T^ax- Parameter k is associated 
with penalty severity: larger values mean stronger penalties. Fitness function is 
adaptive: at the beginning k = 1 and every 10 generations there is a check - if 
more than a percent of individuals are infeasible then k is increased by 0.1. 

4.2 Initialization 

At the beginning initial population of Psize random paths is generated. Con¬ 
struction of each solution begins with adding start vertex s. Afterwards, random 
vertices are iteratively added at the end of the path as long as further insertions 
are possible without violating T^ax constraint (including cost of traveling back 
to end vertex). 
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4.3 Crossover 

During crossover ’’ different pairs of parents are randomly selected from 

the population (cp - crossover probability) and each pair undergoes crossover 
procedure. A specialized 2-point crossover operator was used. At the beginning 
an ordered set S of common vertices for both parents is determined. Vertices or¬ 
der in the set is the same as in first parent. Next, chromosome fragments between 
two successive vertices in S are exchanged and two offspring are created. If any 
duplicates appear in offspring individuals outside of exchanged fragments they 
are immediately removed. Crossover can be performed in two ways depending on 
how exchanged fragments are chosen. Random crossover variant chooses crossing 
points randomly while heuristic crossover maximizes fitness function of the fitter 
offspring (all exchange options are checked). Created offspring can exceed T^ax 
limit. The probability of using heuristic crossover is indicated by parameter Ch 
while the probability of using random version is 1 — c?i. In Figure I there is an 
example of crossover. 


parents 

parent A: (16, 3, 7, 4, 5, 9,12,11,14, 1) 
parent B: (16, 6, 8, 2, 4, 13,10, 9,15, 1) 

crossover variant I crossover variant II 

child 1: (16, 6, 8, 2, 4, 5, 9, 12, 11, 14, 1) child 1: (16, 3, 7, 4, 13, 10, 9, 12, 11, 14, 1) 
Child 2: (16, 3, 7, 4, 13, 10, 9, 15, 1) child 2: (16, 6, 8, 2, 4, 5, 9, 15, 1) 

crossover variant III 
child 1: (16, 3, 7, 4, 5, 9, 15, 1) 
child 2: (16, 6, 8, 2, 4, 13, 10, 9, 12, 11, 14, 1) 

Fig. 1: Example of 2-point crossover. Given two parents which have 4 common vertices 
(16, 4, 9, 1) there are three possible fragment exchanges between successive points. 
Variant I shows children created by fragments exchange between vertices 16 and 4, 
variant II is created by segments exchange between points 4 and 9 while variant III 
shows offspring individuals created by exchange of fragments between vertices 9 and 1. 


4.4 Selection 

In the algorithm no parent selection is used. Instead, survival selection is applied 
in the form of deterministic crowding [3]. After crossover procedure ends each 
offspring competes with one of its parents and the fitter individual is chosen to 
the next generation. To preserve population diversity competition is performed 
between more similar child-parent pairs. Pair arrangement is determined by dis¬ 
tance measure. In this problem two distance measures are applied: edit distance 
(operations insert and delete) and cardinality of symmetric difference of vertex 
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sets. The second is computed initially - if it suggests that child-parent pairs 
should be swapped then more complex procedure of edit distance is calculated. 

4.5 Mutation 

During mutation phase Psize • Wp individuals are randomly selected from the 
population (wp is mutation probability). First, each chosen individual under¬ 
goes 2-opt procedure. This procedure tries to replace two non-adjacent edges 
in a path by two other edges in order to shorten the path as much as possi¬ 
ble. Given a path {vi,V 2 , ■■■,Vi,Vi+i, ...,Vk) 2-opt procedure replaces 

edges (vi, Vi+i) and {vj,Vj+i) by edges {vi, vj) and (vi+i,Vj+i) and reverses path 
fragment between and Vj. The resulting path is: 

(ui, V2 , ..., Vij Vj , Uj _i..., . 

Afterwards one vertex insertion or one vertex deletion is carried out. Proba¬ 
bilities of insertion and deletion are the same (0.5). Both insertion and deletion 
can be either random or heuristic (local search). Heuristic insertion from all 
possible non-included vertices (and all insertion places) chooses the option that 
maximizes ratio of vertex profit to path travel time increase. Random insertion 
chooses new vertex randomly but insertion place is chosen to minimize travel 
time increase. Heuristic deletion chooses the vertex that minimizes ratio of vertex 
profit to path travel time decrease. During mutation T^ax limit can be exceeded. 
The ratio between local search and randomness is determined by parameter mh 
(analogically to Ch in crossover phase). 

4.6 Disturb 

Disturb procedure can be perceived as a different type of mutation which causes 
bigger changes in routes. The operator chooses Psize ' dp individuals from the 
population {dp is disturb probability). Each chosen individual is disturbed - a 
path fragment of random length (bnt containing at most 10 percent of all path 
vertices) is removed. Two types of disturb are possible: first version removes a 
random fragment of random length while heuristic version works analogically 
to deletion mutation - from all possible fragments of a given length it chooses 
the option minimizing ratio of fragment profit to path travel time decrease. 
Specificity of this operator is determined by parameter dh- 

4.7 Adaptation to time dependency 

TDOP algorithm operates on time-dependent networks and for this reason each 
path modification is associated with recalculation of edge travel times after the 
modification point. This can be particularly time-consuming in case of heuristic 
operators. For this reason some modifications were applied to these parts. 
Vertex insertion: After inserting new vertex all edge travel times after inser¬ 
tion place have to be recalculated. Insertion operators of the algorithm choose 
the best insertion place for an included vertex (minimizing travel time increase). 
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Trying each insertion place and recalculating weights is time consuming. How¬ 
ever, it is possible to determine the best insertion place in 0{N) run time (TV - 
number of nodes in a path). The purpose of insertion is to minimize path travel 
time increase but it can be replaced by an equivalent goal of earliest possible 
arrival into the last vertex in a modified path. It can be computed in one loop 
over the path. In the i-th iteration we compare two arrival times into vertex i: 
one resulting from insertion new vertex directly before i and one resulting from 
best insertion place chosen so far. If current insertion place is better than previ¬ 
ous best (arrival time into vertex i is earlier) we update our best choice. 

2-opt: The purpose of 2-opt is to decrease path travel time as much as possi¬ 
ble by edges replacement. Each edge exchange is associated with recalculating 
travel times in a large fragment of a path (part of it is reversed). In order to 
decrease execution time of 2-opt some modifications were applied - instead of 
time-consuming repeated recalculations of travel times the procedure tries to es¬ 
timate total travel time of reversed path fragments. At first the initial travel time 
of the path (before changes) is divided into equal time intervals. Afterwards some 
reversed path fragments (starting at different times) are precomputed. Estima¬ 
tion is fast and base on finding an appropriate precomputed path (explanation 
in Figure 2). Exact calculations are performed only if estimation signals that 
travel time can be shortened - it significantly reduces computations. 



Fig. 2: Estimating travel time of reversed fragment of a path. The path is {vi,V 2 ,..., fio) 
and it starts at time t = 0 and ends at time 80. Time was divided into D — 4 equal 
intervals and all possible vertices in reversed fragments are ug, vg, ... V 2 . The figure can 
be treated as 2-dimensional array and in each cell there is a precomputed path fragment 
visiting a given vertex (determined by column) at a given time interval (determined 
by row). For example in order to estimate travel time of fragment {vr, vg, vs) starting 
at time t = 27 we should find an appropriate column (117) and row (interval 20-40) 
and route fragment associated with it (pointed by arrow). Afterwards we compute the 
difference of arrival times between V 5 and V 7 (deltaT), which is the estimation. 





VIII 


5 Experimental results 

5.1 Parameters 

Experiments were conducted on a computer with Intel Core i7 3.5 GHz processor 
and algorithms were implemented in CH—h. In Table 1 there are values of all 
parameters of proposed evolutionary algorithm. Population size was set to 100 or 
30. The hrst population size {Psize = 100) is derived from original OP algorithm, 
which performed computations on larger networks. The smaller population size 
{Psize = 30) was introduced in order to reduce execution time in smaller TDOP 
instances. Parameters associated with generations were adjusted to population 
size - their values allow algorithm to fully converge. The remaining parameters of 
the algorithm were determined in a process of automatic tuning. The author used 
Parameters ITS [14] meta-algorithm to calibrate his method. The parameters of 
evolutionary algorithm were tuned on classic OP calibration instances but proved 
to be applicable to TDOP instances as well. More details about tuning can be 
found in one of author’s previous papers on the OP [27]. 


Table 1: Parameters of the evolutionary algorithm 


Par am. 

Value 

Description 

Par am. 

Value 

Description 

P . 

Size 

100/30 

population size 

nip 

1 

mutation probab. 

N, 

5000/1500 

max. generations number 

Cp 

1 

crossover probab. 

c. 

500/150 

max. generations number 
without improvement 

dp 

0.01 

disturb probab. 

mu 

1 

mutation heuristic coeff. 

a 

90 

max. percentage of 
infeasible solutions 

Ch 

0.8 

crossover heuristic coeff. 

dh 

0,8 

disturb heuristic coeff. 


5.2 Test instances 

TDOP benchmark instances were introduced by Verbeeck et al [22]. These test 
instances model street traffic. The benchmark authors introduced 7 networks 
containing 21-102 vertices. Maximum travel times {Tmax) were between 5 and 
14 hours and all paths started at 7 in the morning in node 1 and ended in last 
node: N. There are three classes of benchmark networks: 

— Class I contains the smallest networks (21-32 vertices). These are original 
TDOP instances which were solved optimally by benchmark authors. 

— Class II contains all benchmark networks - they were slightly modified so 
that optimal results were the same as in static OP network. Details in [22]. 

— Class III instances were developed from previous classes by applying time 
discretization in order to compare author’s algorithm with Gunawan’s meta¬ 
heuristics [23]. Original start time was 7 am and time limit 14 hours. 
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Table 2: Results for class I instances. N is network size. 


N 

Tmax 

1 EVOlOO 

EVO30 

ACS 

gap 

opt. 

result 

N 

Tmax 

EVOlOO 

EVO30 

ACS 

gap 

opt. 

result 

time 

gap 

time 

gap 

time 

gap 

time 

gap 

32 

5 

0.8 

0 

0.1 

0 

0 

115 

21 

9 

0.6 

0 

0.1 

1.5 

0 

260 

32 

6 

1.0 

0 

0.1 

0 

0 

135 

21 

10 

0.7 

0 

0.1 

0 

0 

310 

32 

7 

1.1 

0 

0.1 

0 

0 

160 

21 

11 

0.7 

0 

0.1 

0 

0 

340 

32 

8 

1.2 

0 

0.1 

0 

2.2 

185 

21 

12 

0.8 

0 

0.1 

0.2 

0 

375 

32 

9 

1.4 

0 

0.1 

0 

0.5 

210 

21 

13 

0.9 

0 

0.1 

0 

0 

425 

32 

10 

1.6 

0 

0.2 

0 

0.4 

230 

33 

5.5 

0.8 

0 

0.1 

0 

4.3 

370 

32 

11 

1.6 

0 

0.2 

0 

0 

250 

33 

6.5 

0.9 

0 

0.1 

0 

0 

420 

32 

12 

1.8 

0 

0.2 

0 

1.9 

270 

33 

7.5 

1.1 

0 

0.1 

0 

4 

500 

21 

5 

0.5 

0 

0.1 

0 

0 

100 

33 

8.5 

1.2 

0 

0.1 

0 

2.9 

560 

21 

6 

0.5 

0 

0.1 

0 

0 

150 

33 

9.5 

1.4 

0 

0.1 

0 

1.9 

620 

21 

7 

0.6 

0 

0.1 

0 

0 

195 

33 

10.5 

1.6 

0 

0.1 

0 

0.3 

650 

21 

8 

0.6 

0 

0.1 

0.3 

0 

220 

33 

11.5 

1.8 

0 

0.2 

0 

1.2 

690 










avg. 

1.0 

0 

0.1 

0.1 

0.8 

322.5 


5.3 Results 

The author tested the evolutionary algorithm (EVO) on these instances and com¬ 
pared it to optimal solutions, ant-colony metaheuristic (ACS) [22] and Adaptive 
ILS [23]. Gaps were expressed in percent and computed as 100 • (1 — where 
res is result obtained by a metaheuristic (average from 30 runs in case of EVO) 
and opt is profit of optimal solution. Execution time is given in seconds. 

In Table 2 results of class I instances are given: performance of two evolu¬ 
tionary algorithm versions {Psize of 100 and 30) is presented and compared to 
ACS and optimal solutions. EVOlOO obtains optimal solutions for all test cases 
and is on average 0.8 percent better than ACS. The biggest obtained gap is over 
4 percent. EVO30 is only 0.1 percent worse than EVOlOO and is much faster. 
Only in one case EVO30 is much worse than EVOlOO (over 1.5%). 

In Table 3 results of class II instances are presented. EVOlOO again achieves 
optimal solutions for most test cases and average gap is close to 0 percent. On 
average it is 1.4 percent better than ACS and for some instances the difference 
is more than 5 percent. EVOlOO clearly outperforms ACS especially in case of 
some larger instances. EVO30 is on average only 0.1 percent worse than EVOlOO 
but for some of largest test cases (100-102 nodes) the difference are larger (0.5- 
1.1 percent). Execution time of proposed methods is also acceptable. 

In Table 4 results for class III instances are compared. The proposed method 
was compared with Gunawan’s Adaptive ILS method (which achieved generally 
the best results among heuristics tested by them). Both versions of EVO are on 
average 1.7-2 percent better than Adaptive ILS and for most instances achieve 
optimal results or close to best known solutions (average gap 0.2-0.5 percent). 
EVO30 performance is generally slightly worse than EVOlOO but it is about 1- 
1.5 percent worse than EVOlOO for a few larger networks. Gunawan’s heuristics 
had execution time limit of 1 second (this limit was also used in case of EVO30). 
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Table 3: Results for class II instances. 


N 

Tmax 

EVOlOO 

EVO30 

ACS 

gap 

opt. 

result 

N 

Tmax 

EVOlOO 

EVO30 

ACS 

gap 

opt. 

result 

time 

gap 

time 

gap 

time 

gap 

time 

gap 

32 

5 

0.9 

0 

0.1 

0 

0 

135 

100 

5 

2.3 

0 

0.3 

0.5 

0.8 

486 

32 

6 

1.0 

0 

0.1 

0 

0 

165 

100 

6 

2.7 

0 

0.3 

0.1 

1.7 

590 

32 

7 

1.3 

0 

0.1 

0 

0 

185 

100 

7 

4.5 

0.2 

0.4 

0.3 

2.5 

679 

32 

8 

1.3 

0 

0.1 

0 

0 

210 

100 

8 

4.0 

0 

0.5 

0.4 

4.5 

771 

32 

9 

1.5 

0 

0.1 

0 

0 

240 

100 

9 

5.0 

0 

0.7 

0.6 

4.2 

853 

32 

10 

1.7 

0 

0.1 

0 

1.5 

260 

100 

10 

6.1 

0.02 

0.8 

1.1 

3 

932 

32 

11 

1.7 

0 

0.2 

0 

0 

275 

100 

11 

6.4 

0 

0.9 

0.2 

0.9 

1007 

32 

12 

1.9 

0 

0.2 

0 

0 

285 

100 

12 

8.8 

0 

1.1 

0.1 

3 

1083 

32 

13 

1.9 

0 

0.2 

0 

0 

285 

100 

13 

7.9 

0 

1.2 

0.1 

1.7 

1147 

21 

5 

0.6 

0 

0.1 

0 

0 

165 

100 

14 

9.9 

0 

1.2 

0.6 

5.3 

1198 

21 

6 

0.7 

0 

0.1 

0 

0 

200 

64 

6.5 

2.6 

0 

0.3 

0 

0.3 

870 

21 

7 

0.6 

0 

0.1 

0 

0 

225 

64 

7 

2.9 

0 

0.4 

0.04 

0.6 

930 

21 

8 

0.7 

0 

0.1 

0 

0 

275 

64 

8 

3.4 

0 

0.4 

0.1 

0.9 

1056 

21 

9 

0.7 

0 

0.1 

0 

0 

315 

64 

9 

3.9 

0 

0.4 

0 

0.8 

1152 

21 

10 

0.8 

0 

0.1 

0 

0 

375 

64 

10 

4.5 

0 

0.5 

0.02 

1.1 

1236 

21 

11 

0.8 

0 

0.1 

0 

0 

415 

64 

11 

6.1 

0 

0.8 

0.4 

1.8 

1308 

21 

12 

0.8 

0 

0.1 

0 

0 

440 

64 

12 

6.2 

0 

0.7 

0.01 

1.5 

1344 

21 

13 

0.9 

0 

0.1 

0 

0 

450 

64 

13 

6.0 

0 

0.6 

0 

0 

1344 

33 

5.5 

0.9 

0 

0.1 

0 

0 

430 

64 

14 

6.1 

0 

0.6 

0 

0 

1344 

33 

6.5 

1.0 

0 

0.1 

0 

0 

490 

102 

4 

2.3 

0 

0.3 

0 

5.3 

532 

33 

7.5 

1.3 

0 

0.1 

0 

0 

550 

102 

5 

2.9 

0 

0.4 

0 

6.1 

648 

33 

8.5 

1.3 

0 

0.1 

0 

0 

590 

102 

6 

3.9 

0 

0.5 

0.2 

1.1 

774 

33 

9.5 

1.5 

0 

0.2 

0 

0 

630 

102 

7 

5.0 

0 

0.6 

0 

1.8 

884 

33 

10.5 

1.7 

0 

0.1 

0 

0 

680 

102 

8 

6.5 

0 

0.8 

0.1 

6.1 

994 

33 

11.5 

1.8 

0 

0.2 

0 

0 

730 

102 

9 

7.2 

0 

1 

0.1 

2.5 

1090 

33 

12.5 

1.9 

0 

0.2 

0 

1 

770 

102 

10 

7.7 

0 

1.1 

0.01 

3.1 

1175 

33 

13.5 

2.0 

0 

0.2 

0 

0.8 

800 

102 

11 

8.6 

0 

1.4 

0 

4.6 

1251 

66 

5.5 

1.5 

0 

0.2 

0 

0.3 

580 

102 

12 

9.6 

0 

1.3 

0 

3.2 

1317 

66 

6 

1.8 

0.03 

0.2 

0.2 

2.3 

650 

102 

13 

15.5 

0.02 

2 

0.5 

5.6 

1368 

66 

7 

2.5 

0.04 

0.3 

0.3 

1.6 

770 


avg. 

3.5 

0.01 

0.4 

0.1 

1.4 

710.7 


6 Conclusions and further research 

The experiments show that proposed evolutionary algorithm with memetic oper¬ 
ators is an effective metaheuristic for the Time-Dependent Orienteering Problem 
(TDOP). The results obtained for TDOP benchmarks by the algorithm are op¬ 
timal or close to optimal in most cases and are better than those obtained by 
other metaheuristics. Author’s method, which was successfully used to tackle 
the classic Orienteering Problem (OP), achieved high quality results for TDOP 
benchmarks as well. Promising results suggest that that proposed method should 
be adapted to other related problems. Therefore further research will concentrate 
on application of the algorithm for other problems from OP family i.e. variants 
with multiple paths (Team OP and Team TDOP) and with time-windows. 
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Table 4: Results for class III instances. Symbol d is discretization step (in minutes). 


d 

N 

EVOlOO 

EVO30 

Ad.ILS 

best 

d 

N 

1 EVOlOO 

EVO30 

Ad.ILS 

best 

time 

gap 

time 

gap 

gap 

result 

time 

gap 

time 

gap 

gap 

result 


32 

1.7 

0 

0.2 

0.0 

0.1 

285 


32 

1.2 

3.8 

0.1 

3.8 

0 

260 


21 

0.8 

0 

0.1 

0 

0 

450 


21 

0.6 

0 

0.1 

0 

2.5 

405 


33 

1.8 

0 

0.2 

0 

0.3 

780 


33 

1.4 

0.3 

0.1 

1.3 

2.5 

680 

1 

100 

16.4 

0.1 

1.0 

0.6 

8.4 

1145 

15 

100 

5.0 

0.1 

0.6 

0.9 

5.9 

864 


66 

6.7 

0 

0.8 

0.2 

1.0 

1513 


66 

1.7 

0 

0.2 

0 

0.6 

870 


64 

7.9 

0 

0.8 

0 

0.5 

1344 


64 

4.9 

0 

0.5 

0.6 

1.6 

1254 


102 

36.6 

0.1 

1.0 

1.4 

6.5 

1372 


102 

9.3 

0 

1.0 

0.5 

4.5 

1068 


32 

1.5 

0 

0.2 

0 

0 

280 


32 

0.9 

0 

0.1 

0.2 

0 

220 


21 

0.7 

0 

0.1 

0 

2.2 

450 


21 

0.5 

0 

0.0 

0 

0 

355 


33 

1.5 

0 

0.1 

0 

2.1 

760 


33 

1.1 

0 

0.1 

0 

0 

590 

5 

100 

9.3 

0 

1.0 

0.6 

7.6 

1032 

30 

100 

2.4 

1.2 

0.3 

1.7 

3.5 

623 


66 

3.2 

0 

0.4 

0.3 

1.0 

1260 


66 

1.6 

0 

0.2 

0 

0.2 

850 


64 

6.1 

0 

0.6 

0.5 

2.9 

1274 


64 

1.5 

0 

0.1 

0 

0 

768 


102 

14.4 

0 

1.0 

1.5 

6.8 

1244 


102 

2.6 

0.1 

0.3 

0.3 

1.1 

690 










avg. 

5.1 

0.2 

0.4 

0.5 

2.2 

810.2 
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